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The size of the emission volume of pions and kaons, r wn and r^x, were measured 
in the hadronic Z° decays through the two-pion and two-kaon Bose-Einstein correlations 
near threshold. Recently the emitter size of the two identical baryons, AA(AA), was 
evaluated from the dependence of the fraction of the spin S=l state on the energy near 
threshold where it is affected by the Pauli exclusion principle. Here we show that the r 
dependence on the particle masses, namely the hierarchy r wn > tkk > r AA observed in 
the Z° hadronic decays, is well described in terms of the Heisenberg uncertainty principle. 
A good description can also be obtained via the virial theorem when applied to a general 
QCD potential. Other available approaches are also discussed. 
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1 Introduction 



It has been known for more than four decades that the correlation function of two identical par- 
ticles allows us to learn about the size of the particle emission region (the Hanbury Brown and 
Twiss effect [0). The interference effect originated from Bose-Einstein (Fermi-Dirac) statistics 
leads to an enhancement (reduction) of the two-particle correlation function which occurs when 
two identical particles are bosons (fermions) and have almost equal momenta. The relation to 
the size of the emission region is often given by the well known formula for the correlation 
function of two identical particles with four momenta g« (i = 1, 2) and Q = J—(qi — (h) 2 -, 
namely 

atot ^k S ° {Q) = 1 ± |jR(g)|2 ' (1) 
where the sign + (— ) corresponds to bosons (fermions) and R(Q) is the normalised Fourier 
transform of p, the source density || 

iQ) = Wpjx) ' { ) 

This quantity is often parametrised in terms of a "source radius" r, related to the size of the 
emission region, and a "chaoticity parameter " A, which measures the strength of the effect. 
From Eq. (§) one derives that: 



\R(Q)\ 2 = Xe~ r2Q2 , (3) 

assuming a spherical emission volume with a Gaussian density distribution. The study of the 
Bose-Einstein Correlations (BEC) of identical boson pairs, is carried out by forming the ratio 
of the data correlation function to the same correlation of data or Monte Carlo (MC) sample 
having all the correlations, like those arising from resonances apart from those due to BEC. Since 
the original study of Goldhaber et al. || in 1959 and 1960, the BEC of pion pairs were studied 
in a large variety of interactions over a wide range of energies Q. In addition, studies of BEC of 
kaon pairs have also been reported and values for their emitter dimension have been evaluated. 
The study of this dimension has recently been extended to pair of baryons, namely the AA(AA) 
systems using the Pauli exclusion principle. From these experimental results an interesting 
hierarchy of the size of the emission volumes is observed, namely r n7T > tkk > ?"aa- This 
volume size dependence on the hadron mass may give a new insight to the non-perturbative 
QCD and the hadronisation process. In this work we investigate this dependence using ttti, 
KK and AA pair samples present in the multi-hadron Z° decays. In Section ^| we summarise 
the current knowledge of the BEC parameters obtained from the studies of the hadronic Z° 
final states and describe our procedure to extract representative taa values. In Section [| we 
show that by using the Heisenberg uncertainty principle one is able to describe successfully the 
r dependence on the particle mass value. In Section f| we extract the mass dependence of the 
size of the interaction volume which follows from the virial theorem by using a general QCD 
potential. Finally the information on the interaction character that can be derived from the 
identical particle correlations is summarised and discussed in Section 5. 
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2 Extraction of the two-particle emitter size 



During the last decade values for r nn have been reported by the ALEPH ||, DELPHI [[J and 
OPAL |7] collaborations from measurements in the hadronic Z° decay events. These r values 
are summarised in Table To notice is the fact that the r values are rather sensitive to the 
choice of the reference sample which was either the data 7r + 7r~ system or a sample of two-pions 
originating from two different events, the so called mixed event sample. The ALEPH and 
DELPHI collaborations have therefore opted to quote their final result as the average of the 
values obtained from these two methods. OPAL chose as the best r value the one obtained 
with a 7r + 7r~ data reference sample. The result obtained with a MC generated reference sample 
was used for the estimation of the systematic error. To be on equal footing with the ALEPH 
and DELPHI collaboration, we also present in Table 1 the average OPAL value obtained from 
their two methods. Finally we also list in the last row of Table |l| the over-all averaged r V7T value 
of 0.74 ± 0.01 ± 0.14 fm where the first error and the second one are respectively the statistical 
and systematic uncertainties. 



7W (fm) 

Experiment Method I Method II I & II Combined 
ALEPH § 0.80 ±0.04 0.50 ± 0.02 0.65 ± 0.04 ± 0.16 
DELPHI [| 0.82 ±0.03 0.42 ± 0.04 0.62 ± 0.04 ± 0.18 
OPAL 0.96 ±0.01 0.79 ±0.02 0.88 ± 0.02 ± 0.09 
LEP Average 0.90 ± 0.01 0.62 ± 0.01 0.74 ± 0.01 ± 0.14 



Table 1: The two-pion emitter size obtained from BEC studies using the hadronic Z° decays 
carried out with two choices for the reference sample. In method I the reference sample was the 
data 7r + 7r~ system. In method II the reference sample was created by the event mixing technique 
or Monte Carlo generated events. In the last row the weighted average values are presented. 
The first and the second errors are respectively the statistical and systematic uncertainties. 
The systematic error is taken as half the difference between the values of the two methods. 



Several BEC studies have been reported for the pairs of kaons. The BEC parameters obtained 
for the KgKg system do suffer from the fact that part of the observed low Q enhancement may 
originate from the decay of the / (980) — > K Q ~K° . On the other hand the BEC of the K ± K ± 
pairs are free of this problem. In the following whenever we refer to a kaon pair we mean the 
K K systems. Using the Z° hadronic decays, DELPHI || has measured the BEC of these 
pairs and obtained results for the emitter dimension tkk and the strength Xkk- 

r KK = 0.48 ± 0M(stat.) ± 0.07{syst.) fm and \ KK = 0.82 ± 0.11(stot.) ± 0.25{syst.) 

A method for the measurement of the emitter dimension of two identical baryon systems in 
multi-hadron final state has recently been proposed by G. Alexander and H.J. Lipkin ||. This 
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method, which is applicable to pairs of identical spin 1/2 baryons which decay weakly, uses 
the Fermi-Dirac statistics properties which suppress the over-all spin S=l state of the identical 
di-baryon system near their threshold (Pauli exclusion principle). 

To apply e.g. this method to AA or AA pairs, the hyperon-pair is first transformed to their 
common centre of mass (CM) system and then each A decay proton is transformed to its 
parent hyperon CM system. The distribution of the cosine angle between these two protons, 
here denoted by y*, depends on the fraction of the S=l (or S=0) in the data. Specifically the 
following y* distributions are expected for pure S=0 and S=l states of the A pair 

dN/dy*\ s=0 = 1 + {-I) 13 ' 2 ■ a 2 A ■ y* and dN/dy*\ s=1 = 1 - {-l) B ' 2 ■ a\ ■ y*/3 , 



where a A =0.642±0.013 |TIJ is the A — ► pn~ decay asymmetry parameter arising from parity 
violation and B is the baryon number of the di-hyperon system. These distributions are in- 
dependent of the orbital angular momentum and are valid as long the A's are non-relativistic 
in their CM system. Therefore in this method the spin composition analyses are restricted to 
relatively small Q values^. Defining e as the fraction of the S=l contribution to the di-hyperon 
pairs, it is equal to 

(S = 1) 
(S = 1) + (S = 0) ' 

so that the following function: 



dN/dy* = f BG + (1 - f BG ) • {(1 - e) ■ dN/dy*\ s=0 + e ■ dN/dy*\ s=1 } , 

can be fitted to the data at every Q values. Here Jbg is the background fraction in the data. 
As in BEC analyses the possible, if at all, final states interactions are neglected. 

The fraction of the S=l spin content, e, of the AA(AA) system has been measured as a func- 



tion of Q in ~ 4x 10 hadronic Z decays per experiment at LEP by the OPAL ||11|| , ALEPH [[12 
and DELPHI [13j collaborations. A decrease in e has been observed as Q approached zero in 
contrast to the results obtained for the AA where e was found to be equal to 0.75 down to very 
low Q values. The taa values deduced by the OPAL and DELPHI collaborations are : 



opal m 



r AA = 0.19 



-0.37 
-0.07 



± 0.02 fm 



DELPHI [13] r AA = 0.11 



+0.05 
-0.03 



± 0.01 fm . 



The ALEPH results for the AA spin composition are similar to those of OPAL and DELPHI 
but no attempt has been made to extract from them an r AA value. 

The measured e(Q) values of the three experiments are compiled in Fig. |I|a where the error 
bars are strongly dominated by the statistical uncertainties. The solid line in the same figure 
is the outcome of our unbinned maximum likelihood fit of the function: 

e(Q) = 0.75[1 - e" r iA Q2 ] , (4) 
4 Assuming that no exotic AA resonances exist in this Q region. 
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Figure 1: (a) The S = 1 fraction, e, of the AA(AA) pairs measured as a function of Q by 
the ALEPH [|T2J, DELPHI (I3j and OPAL |TT| collaborations. The solid line represents the fit 
results of Eq. (f|) to the data points, (b) The A% 2 = x 2 — Xmin dependence on taa as obtained 
from the fit to the combined LEP data. 



to the data points plotted in the Fig la. The factor 0.75 appearing in Eq. (|4|) represents 
the statistical spin mixture which is proportional to 2S + 1. From this over- all fit, with a 
X 2 /d.o.f. = 0.4, the dimension of the AA emitter and its errors are found to be 

r AA = 0.14+°^ f m . 

The Ax 2 = X 2 ~ Xmin behaviour of this fit is shown in Fig. |I|b. The ALEPH collaboration has 
also used an alternative method to measure r A A- In that method one constructs, similar to the 
BEC studies, a correlation function of the type 

C^{Q) = NiUQ)/N r \ A f (Q) , 

where the numerator is the data distribution and the denominator is the distribution of a 
reference sample void of the Fermi-Dirac statistics. This C^ A (Q) correlation is expected to 
decrease at low Q values due to the onset of the Pauli exclusion principle. From this analysis 
ALEPH evaluates taa to be (0.09 — 0.10)±0.02 fm in perfect agreement with the values obtained 
from the spin composition analyses. 

In Fig. § we present the experimental results for r obtained in the analyses of the hadronic 
Z° decays. The DELPHI values, which is the only LEP experiment that measured Tkk-, are 
shown by circles. The LEP averages of the measured r n7T and r A A values are also shown by 
triangle symbols and are seen to be in good agreement with DELPHI ones. The large error 
associated with r nn reflects the systematic uncertainty due to the strong dependence of the 
experimental results on the choice of the reference sample. In contrast to this situation, a 
relative small systematic error is associated with taa since in this case the spin analyses do not 
rely on a reference sample. In spite of the relative large errors a general non-trivial trend can 
be observed in the r dependence on the particle mass, namely that: 

' 7T7T > r KK > r AA . 
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Hadron mass (GeV) 



Figure 2: The emitter radius r as a function of the hadron mass m. The triangles are the 
DELPHI results and the circles are the averages of the measured values at LEP. For clarity the 
points are plotted at slightly displaced mass values. The error bars correspond to the statistical 
and systematic errors added in quadrature. The thin solid line represents the expectation from 
the Heisenberg uncertainty relations setting At = 10~ 24 seconds. The upper and the lower 
dashed lines correspond to At = 1.5 x 10 -24 and 0.5 x 1CT 24 seconds respectively. The thick 
solid line represents the dependence of r on m as expected from a QCD potential given by Eq. 

m 

3 The emission volume size and the uncertainty relations 



The maximum of the BEC enhancement of two identical bosons of mass m occurs when Q — > 
which also means that the three vector momentum difference of the bosons approaches zero. 
This motivated the interest to link the BEC effect to the uncertainty principle |[14|| . From the 
Heisenberg uncertainty relations we have that 



ApAr = 2 ji vr = m vr = he , (5) 

where m and v are the hadron mass and its velocity. Here \x is the reduced mass of the di- 
hadron system and r is the distance between them. In Eq. (|5]) the momentum Ap is measured 
in GeV, Ar = r is given in fermi units and he = 0.197 GeV fm. From this, one obtains 



he he 

r = — = — . (6) 

mv p 

Simultaneously we also use the uncertainty relation expressed in terms of time and energy 
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v 2 

AEAt = f- At = K , (7) 

m 



where the energy is given in GeV and At in seconds. Thus we have 



p —hm/At so that p = yhm/At. (8) 
Inserting this expression for p in Eq. ([]) we finally obtain 



hc/Jh/At c ^nKt , n . 

rim) = —= — = — ■=— . (9) 

In the following we take At = 10 -24 seconds as a representative time scale of the strong 
interaction which is assumed to be independent of the hadron identity and its mass. Thus 
r(m) = A/yftK with A = 0.243 fm GeV 1 / 2 . 

In Fig. |2| we show by the thin solid line the r dependence on the mass as calculated by Eq. 
(|9|) setting At = 10 -24 seconds. The sensitivity of r{m) on the value of At is illustrated by 
the dashed lines in the figure. The upper and the lower dashed lines correspond respectively 
to the values At = 1.5 x 1CT 24 and 0.5 x 10~ 24 seconds. A fit of Eq. (g) to the three data 
points plotted in Fig. || yields At = (1.2 ± 0.3) x 10" 24 seconds with a \ 2 /d.o.f. = 4.2 which 
corresponds to a probability of ~12%. As seen from Fig. |2], the dependence of r on m as given 
in Eq. (|9]) describes, within errors, very well the mass dependence of the measured emitter size 
and offers a simple explanation for the measured hierarchy r n7T > r^K > r AA- 



4 The emission volume size and the virial theorem 



The previous section shows that the observed hierarchy r nw > r KK > taa has a natural and 
a rather general explanation. Here we would like to demonstrate that such a dependence could 
shed light on the character of the "soft" interaction or, in other words, on the non-perturbative 
QCD which is responsible for the interaction at long distances. Let us continue to use the 
semi-classical approximation as in section ^, which means that the angular momentum 

£ = \pi - p 2 \b t = 2pb t w he, (10) 



where bt is the impact parameter. To estimate the value of p we use the virial theorem |1| 



which yields a general connection between the average values of the kinetic and the potential 
energies, namely 

2(T t ) = (b t -V t V(r)), (11) 

where t denotes the transverse direction. Since the motion in the transverse direction is always 
finite, we can safely use the virial theorem. Substituting in Eq. ([□]) the kinetic energy by its 
relation to the momentum, T t = p 2 /m, we obtain 

(p 2 t ) = m(f t ■ V t V(r) ) , (12) 



6 



where we denote by r the distance between two particles with r t being equal to 2b t . 

From Eq. ([12]) one can see that the relation between {p 2 ) and r depends crucially on the r 
dependence of the potential energy V(r). Here we take V(r) to be independent of the interacting 
particles (quarks) mass. This assumption is also in accordance with the Local Parton Hadron 
Duality [16[ concept, which states that one can consider the production of hadron as interaction 
of partons, i.e. gluons and quarks, in spite of the unknown mechanism of hadronisation. This 
approach is supported experimentally by the single and double inclusive hadron production in 
e + e~ annihilation |16| . Since the interaction between partons is independent of the produced 
hadron mass it follows that dV(r)/dm = 0. Based on this assumption and using Eq. (|It]) and 
Eq. (0), we derive an equation for the typical size of the emission volume, namely 

(fir) 2 

r 2 ( r • VV(r)) « . (13) 

m 

To illustrate the kind of information that can be extracted from Eq. (|T3"|), we evaluate the 
typical size of the particle emission region from two very different potential energy V(r) forms: 



1. V(r) = N c a s hc/ (tt r) which models the QCD induced interaction at short distances. In 
this case, from Eq. (|13"|) one obtains r oc he/ {as fn). 

2. V(r) = Kr which corresponds to the confinement potential. For this potential Eq. ( THf ) 
leads to r oc l/ms. 



> 
O 




r(fm) 



Figure 3: The behaviour of three factors, r ■ VV, r 2 (f ■ V) and r 2 which contribute to the 
dependence of 1/m on the radius r when using the potential V(r) defined in Eq. (0). 



It is interesting to note that the experimental data plotted in Fig. 2 are well described by 
Eq. (^) when applied to the general potential form 

4 otsfic 
3 r ' 



V(r) 



K r 



(14) 
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which is widely used to derive the wave functions and decay constants of hadrons [17]. This 
can be seen by comparing in Fig. 2 the thick solid line with the experimental data on the radii 
derived from the identical particles correlations. This line was obtained by setting Eq. ( |i~3"D to 
an equality and evaluating Eq. (|T4] ) with the V(r) parameter values of k — 0.14 GeV 2 = 0.70 
GeV/fm and a s = 2tc/9 \n(5 + 7/r) with 5 = 2 and 7 = 1.87 GeV -1 = 0.37 fm obtained from 



the hadron wave functions and decay constants |T7 . 

Here it should be stressed that the simple r oc l/^/m dependence obtained from the 
uncertainty principle does not contradict the virial theorem estimate. The fact that for the 
QCD potential defined by Eq. (|T4"D, the factor (r ■ V t V{r) ) is essentially constant (see Fig. 
3) over the r range 0.14 to 1.0 fm, means again that r oc 1/ ' \fm. The proportionality factor 
then sets the scale of distances in the non-perturbative QCD. 



5 Summary and Discussion 

The study reported here was motivated by the recent experimental results on the emission 
volume size derived from the identical particle correlations measurements. The information 
that this size is seen to decrease as the mass of the hadron increases, is shown here to be very 
valuable for the understanding of the non-perturbative QCD mechanism. 

A good description of the r dependence on the hadron mass can already be obtained from the 
uncertainty principle. However this requires the acceptance of two, not trivial and not obvious 
assumptions, namely: 



1. At does not depend on the mass value of the produced particles. As is shown below, this 
assumption is inconsistent with one of the models for the non-perturbative description of 
the confinement forces, namely, with the string approach. 

2. AE depends only on the kinetic energy of the produced particle. It means that the 
interaction (potential energy) is rather small. 

3. Ar ~ r which is in the spirit with the uncertainty principle used here. 



The second approach followed here is based on the virial theorem and uses the Local Parton 
Hadron Duality (T^[, which is the only argument why we can consider the potential energy 
to be independent of the mass of the produced particles. It looks rather impressive that the 
phenomenological potential, whose parameters were taken from the characteristics of the hadron 
wave functions and decay constants, is able to describe so successfully the experimental data. 
This gives us therefore the hope that the experimental data on the radius of the emission volume 
would be useful in any discussion of the theoretical and/or phenomenological description of the 
non-perturbative QCD forces responsible for confinement of quarks and gluons. 

The decrease of the emitter size with the mass of the produced hadrons limits the models 
which can be applied to the non-perturbative QCD. For example, consider the string approach 
to the non-perturbative QCD, which is also the base of the Lund Model of hadronisation 
described in detail in f|Tj|. In this approach the formation time which we need to create a 
particle with mass m, as well as the corresponding distance, are proportional to the mass of 
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the produced particle, namely, t(r) oc tu/k where k is the string tension which appears in 
Eq. (|14]) . The experimental data in Fig. 2 shows in fact the opposite trend. Therefore, the 
string based models may have difficulties in accommodating the data findings. On the other 
hand, the perturbative QCD cascade plus the Local Parton Hadron Duality Jl6| approach leads 
tor oc 1/m which reproduces qualitatively the experimental data. These two examples clearly 
illustrate how these measurements are able to select or reject theoretical proposals. For this 
reason more accurate measurements of r W7r , r KK and should be very valuable. 

We would like to stress that the observed hierarchy of the sizes of the emission volume is 
providing us an additional information on the dynamics of interaction to that given by the well 
established experimentally dependence of the effective source sizes on the produced particles 
transverse mass, m t . Experiments with heavy-ions and hadron-hadron reactions as well as 



LEP studies |1| show that r ~ l/^/m~t. In these experiments only Bose-Einstein correlations 
between identical pions have been measured, where r is considered to be proportional to the 



production time ||20|| . It should be stressed that this production time is defined as the invariant 
time elapsed between the creation of a qq pair and the formation of a pion. For large pt an 
hypothesis has been suggested that the production time is larger than the formation time. This 
hypothesis is able to reproduce the experimental dependence of r on m t . It is reasonable to 
assume that the dependence of the size of the source on the mass of the produced particle is 
related to the formation time rather than to the production one. This is a reason why such 
experiments as discussed here can give a new insight to the nonperturbative dynamics of hadron 
production. 
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